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ABSTRACT 

We present optical and near-infrared photometry of the afterglow of the long 
Gamma- Ray Burst GRB 070802 at redshift 2.45 obtained with the ESO/MPI 
2.2 m telescope equipped with the multi-channel imager GROND. Follow-up 
observations in g'r'i' z' and JHKs bands started ~17 min and extended up to 
28 h post burst. We find an increase in brightness of the afterglow at early times, 
which can be explained by the superposition of reverse and forward shock (FS) 
emission or the onset of the afterglow FS. Additionally, we detect a strong broad- 
band absorption feature in the i' band, which we interpret as extinction from the 
redshifted 2175 A bump in the GRB host galaxy. This is one of the first and 
clearest detections of the 2175 A feature at high redshift. It is strong evidence 
for a carbon rich environment, indicating that Milky Way or Large Magellanic 
Cloud like dust was already formed in substantial amounts in a galaxy at z=2.45. 

Subject headings: gamma-rays: bursts — X-rays: individual (GRB 070802) — 
ISM: dust, extinction 



Introduction 



Gamma-Ray Bursts (GRBs) are intense extrag alactic flashes of 7-rays with du rations 



between several tenths to hundreds of seconds (e.g. 



Fishman 


1994; 


Meszaros 


2006) 
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cluster in two di fferent categories o 



long soft bursts (IKouveliotou et al 
radio to X-rays (e.g. 



duration and spectral hardness: short hard bursts versus 



19931). GRBs are fo l lowed by loii ger lasting afterglows in 



Katall994l : iMeszaros fc Reeslll997l : IZhangll2007l ). The enormous energy 



release in 7-rays is unaffected by dust or gas absorption, so GR Bs are detectable out to very 



high redshifts (e.g. iLamb fc ReichartI I2OOOI : iKawai et al.l 120061 ). These characteristics make 



( Savaglio 


2006; 


Fvnbo et al. 


2006; 


Berger et al. 


2006) 



The GRB-SN connection (e.g. 



Zeh et all 12004 : IWooslev fc BloomI l2006|) strongly sup- 



ports the association of long GRBs with the core-collapse of ver y massive stars, a nd therefore 
with regions of high- mass star formation in the host galaxies (jPaczynskil Il998l ) . High-mass 
star forming regions show signatures of signi ficant dust and gas absorption i r i ultra-violet 



and o ptical spectra of bright GRB afterglows (jSavaglioll2006l : iBerger et al.ll2006l : iFynbo et al. 



20071 ). In fact, a large number of bursts do not show optica l afterglows des pite rapid and 
deep optical follow-up observations f e.g. IStratta et al.l 2004 iRol et al.l l2005l ). Natural ex- 
planations f or 'optically dark bursts' (iGroot et al.lll998l ) include dust ex tinction in the hos t 



galaxy (e.g. IPvnbo et allboOll : Ixiose et ahlbooal ) and high redshift fe. g. iGroot et al 



both efficiently suppressing any fiux in the observed optical bands (iRoming et al 
Several bursts with moderate dust absorption have been detected (e.g. 



Kann et al. 



which may represent only the lower end of the host extinction distribution (jSchady et al. 



1998h 



20061) 



20061 



2OO7I ). The detection of significant reddening in a host is strongly instrumentally biased due 
to the lack of rapid near-infrared follow-up observations with large aperture telescopes. In 
particular, one of the biggest shortcomings of all multi-color monitoring of afterglows so far 
was the lack of a simultaneous coverage of the optical to near-infrared (NIR) bands. 

The Gamma-Ray Burst Optical and Near Infrared Detector (GROND) is a 7- channel 
imager prim arily designed for fast follow-up observations of GRB afterglows (iGreiner et al. 



2007bl . I2OO8I ). It enables the detection and identification of GRB afterglows in a broad 
wavelength range (400 nm - 2310 nm). Due to the use of dichroic beamsplitters it is capable 
of simultaneous imaging in seven bands, g'r'i'z' (similar to the Sloan system) and JHKs- 
GROND is mounted on the 2.2 m ESO/MPI telescope on LaSilla/Chile since April 2007. 
The instrument is operated robotically and capable of monitoring the light curve of the 
transient starting from a minimum of a few minutes after burst alert. The fi eld of view of 



the in strument is 10' x 10' in the NIR and 5.4' x 5.4' for the optical bands (IGreiner et al. 



20081 ) . The fast response, medium sized telescope aperture, NIR capabilities and unique 
optical design makes GROND an ideal tool for follow-up observations of GRBs. 

Here we report on first GROND follow-up observations of GRB 070802 (section [2]) and 
derive constraints on the GRB ejecta and its circumburst properties (section [3]) . 
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Observations 



2.1. Swift observations 



The BAT (Burst Alert T elescope) instrument (jBarthelmy et al.ll2005l ) onboard the Swift 
satellite JCehrels et alibood ) triggere d on the long-soft GR B 070802 at To=07:07:25 UTC 
and immediately slewed to the burst (jBarthelmy et al.l 120071 ) . The BAT light curve shows a 
single peak starting at Tn+5 s and ending at T n+50s. There is evidence at the 3a level for 
a precursor at Tq-ISO s JCummings et al.ll2007h . The Tgo for GRB 070802 is 16.1±1.0 s and 
the fiuence in the 15 keV to 150 keV band is 2.8±0.5 x 10 



20071 ). 



erg cm ^ (iCummings et al. 



The X-Ray Telescope (XRT, iBurrows et al.l l2005l ) began follow-up observations of the 
burst field 138 s after the trigger and detected an uncatalogued fading X-ray source at a 
position of Ra(J2000.0)= 02^ 27"^^" 35^76, Dec l(J2000.0)=-55° 31' 38'.'4 with a refined 90% 



coincidence error of 2'.'1 (IMangano et al.l 120071 ) . The XRT light curve decays with a slope 
of a ~ —2 until To-|-500 s and then remains fiat until ~4 ks. Afterwards the light curve 
decays with a power law index of ~ — 1 until it fades below the XRT sensitivity limits 
(jMangano et al. 2007). No bright flares are detected in the XRT light curve. 



The third instrument onboard Swift, the Ultra- Violet Optical Telescope (UVOT. lRoming et al 



20051 ) started observations at Tq+IOO s and did not find any transient sources inside the 
XRT error circle do wn to 19.5 (u filter), 21.3 (b filter) and 21.2 (uvw2 filter) magnitudes 
Jimmler et al.ll2007l ). 



2.2. GROND optical and near infrared observations 

GROND responded to the Swift GRB alert and initiated automated observations on 
the 2""^ of August 2007 at 07:24:09 UTC, starting 7 min 16 s after the Swift trigger and 
16 min 44 s after the onset of the burst. A predefined sequence of exposures with succes- 
sively increasing exposure times were executed, acquiring images in all seven photometric 
bands simultaneously. The observations continued for two nights, after which the afterglow 
had faded below the GROND sensitivity limits. A variable point source in the NIR bands 



(iGreiner et al.ll2007al ) inside the Swift XRT error circle was identified by the GROND data 
reduction pipeline (A. Yolda§ et al. 2008, in preparation). The transient is shown in Fig. [T] 
and its absolute position is calculated to be Ra(J2000.0)=02'' 27^" 35^.68, Decl(2000.0)=- 
55° 31' 38'.'9 with an uncertainty of 0'.'3 compared to 2MASS reference field stars. The 
afterglow was also observed and detected by the Magellan telescope at Las Campanas Ob- 
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servatory (LCO. iBerger &: Murphyl 120071 ) and the Very Large Te lescope (VLT), the latter 
yielding a spectroscopic redshift of z=2.45 (jProchaska et al.l 120071 ). 



In total 2036 NIR images with an integration time of 10 s each, and 56 CCD optical 
frames were obtained with GROND. The CCD integration times scaled with the brightness 
of the transient from 45 s at early times to 10 min when the source had faded. Sky conditions 
were clear with a mean seeing around 1'.'2. All GROND data were obtained at airmasses 
between 1.35 and 1.12. 

Optical and n ear-infrared image reduction and photometry was performed using stan- 
dard IRAF tasks (ITodyl Il993l ). A general model for the point-spread function (PSF) of 
each image was constructed using bright field stars and fitted to the afterglow. Addition- 
ally, aperture photometry was carried out and the results were consistent with the reported 
PSF photometry. Photometric calibration was performed relative to secondary standards 
in the GRB field. During photometric conditions, two s pectrophotornetric standard stars, 
SAl 14-750 and SAl 14-656, both primary Sloan standards (ISmith et al.ll2002l ). were observed 
with GROND. Observations of the GRB field followed within 4 minutes. The magnitudes of 
SA114-750 and SA114-656 wer e transformed to the GROND filter system using their spectra 
and the GROND filter curves (iGreiner et al.ll2008l ). The obtained zeropoints were corrected 
for atmospheric extinction differences and used to calibrate stars in the GRB field, shown 
in Fig. [H The apparent magnitudes of the afterglow were measured with respect to the 
secondary standards reported in Table [H 

Vega magnitudes have been transformed to the AB system using transformation factors 
for the GROND filter system as 6g'=0.01 mag, 6r'=0.15 mag, 6r'=0.39 mag, 6z'=0.52 mag, 
6J=0.91 mag, 6H=1.38 mag and dKs=l.SO mag. All reported afterglow magnitud es are 
corrected for Galactic foreground reddening (i?B_i/=0.026 mag, ISchlegel et al.lll998l ). As- 
suming Rv=3.1 for the Milky Way, this leads to A^=0.01, Ah=0.02, Aj=0.03, A,=0M, 
A=0.05, v4^=0.07 and ^3=0. 10 for the GROND fiher bands. 

The afterglow is detected in all seven GROND bands, however it was too dim to con- 
struct a light curve with reasonable time resolution in the filter bands g' and i'. The light 
curves obtained in the r'z'JHKs bands are presented in Fig. [21 The light curve behaviour is 
dominated by an early rise in brightness, after which it declines with a bump superimposed 
onto the overall decay. The observed variations occur in all five GROND bands, and the 
generic light curve shape is achromatic within the measurement uncertainties. 
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3.1. 



3. Analysis 
The early light curve of GRB 070802 



Wiiers et al. 


1997; 


Sari & Piran 


1999; 


Piran 


2005 



Meszaroa l2006l ) . GRBs produce their prompt emission in 7-rays from internal shocks of an 
ultra-relativistic outflow from a compact source and long-wavelength afterglows from the 
interaction of the ejecta with the circumburst medium. After the prompt internal shock 
phase, the optical afterglow light curve is composed of a superposition of two different emis- 
sion components. The reverse shock (RS) prop agating into the eje cta and the forward shock 



(FS) travelling into the surrounding medium (IZhang et al.ll2003l ). Rapid optical observa- 



tions of the early tran sition phase between prompt and afterglow emission can constrain the 



nature of the outflow (INakar fc Pirarul2004l ). While baryonic ejecta are expected to produce 
an optical flash, that can be associated with a RS, a Poynting flux dominated afterglow 
should preferentially show the FS emission. The hydrodynamical calculations from the flre- 
ball model have succeeded in describing the generic afterglow of GRBs from several minutes 
to days post burst. However, the majority of burs ts do not show brigh t optical flashes and 



apparently lack a strong RS emission component (IRoming et al.l 120061 ). This fact provides 
observational support that the strength of the optical emission from the RS may be weaker 
than previously calculated jBeloborodovlboOsh . 



A baryonic shell expanding into a homogeneous medium is fully described by the shell 
iso tropic equivalent ener gy E, its width A, the initial Lorentz factor Fq and the ISM density 
n (INakar fc Pirarul2004l ). In the thick shell case, the initial Loren tz factor is large, and th e 
swept up circumburst medium decelerates the ejecta efliciently (IKobayashi fc Saril I2OOOI ). 
Here, the RS becomes relativistic in the rest frame of the unshocked ejecta. The critical 
Lorentz factor Fc discriminatin g between relativis tic thick (Fq > Fc) and Newtonian thin 
shell case (Fq < Fc) is given by ( Zhang et al.lboosi ): 



-i/sj^-a/s 



l + z 



3/8 



where £"52 is the isotropic energy equivalent in units of 10^^ erg and T2 is the burst duration 
in units of 100 s. For the thin shell case the Lorentz factor at the crossing time of reverse 
and forward shock is F^; = Fq. By measuring the peak of the RS, one can directly derive Fq 
for the thin shell case: 

F. = V,{T/t,f'^ (2) 
with only a weak dependence on E^2/n. 

The light curve shown in Fig. [2] displays a behaviour consistent with the upper theo- 
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retical predictions. It rises achromatically at early times and reaches a first peak at ~2 ks. 
Afterwards the hght curve dechnes with an power law index of ~-2 until ~3 ks. Between 
3 ks and 10 ks post burst our data are not conclusive. There may be a plateau phase, al- 
though a power law dechne with the late time index cannot be ruled out. At late epochs the 
afterglow follows the canonical power law decay. The data are compatible with a constant 
power law spectrum at all times, although there is an indication for chromatic changes at 
peak brightness. There is no evidence for a jet break in the GROND data out to 90 ks post 
burst. 

We used two different approaches to analyze the light curve quantitatively. Firstly, a 
physical model combining the forward and reverse shock component (RS-FS model), and 
secondly a forward shock model alone (FS model). 

Within the framework of the combined RS-FS model, the first peak can be interpreted 
as the peak of the reverse shock, whereas the possible rebrightening phase around 4 ks is 
relat ed to the FS peak. A complete parametrizat ion of the RS is given by a broken power 
law (IBeuermann et al.lll999uNakar fc Piranll2004l ): 



Kit) 



(3) 



with power law indices of rise (a^ >0) and decline (0^2 ~ 2), a peak time tp^ak = tiic^i/ — 
^r^i/(s' (a'^-aj))^ normalization Fq, and the sharpness of the break s^. 

The FS parametrization used for the light curve analysis is similar to the RS: 

f f f f-i -i/s^ 



-^0 



+ 



t 



(4) 



According to theoretical calculations (ISari et al.l Il998l ). a power law index of the FS with 



a{ ~ 1/2 is expected for the forward shock rise, followed by a shallow decline ~ 3/4— 3/4p 
Here p is the power law index of the energy distribution of the shocked electrons (jSari et al. 



1998 



van Paradijs et al.ll2000l ). The power law indices for both forward and reverse rise and 



decline as well as the peak time have been fitted using the complete multi-color data set. 
The fiux normalizations result from single band fits. The result of the combined fit of the 
RS and FS are consistent with the expected power law decline of the RS with a power law 
index ~ —2 (Tables [2] and [3]). 

To constrain the RS/FS fit, the sharpness parameter of the power law transition (s^ ) and 
the time of the putative FS (t2) had to be fixed to 2.5 and 4.5 ks, respectively. However, this 
affects mainly the late time properties of the afterglow, whereas our analysis is concentrated 
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on the early epoch dominated by the rise of the hght curve. Specifically, the peak of the 
light curve is unaffected by the fixed parameters. 

We also fit the data with the FS broken power law alone. In this case the early rise of 
the afterglow is interpreted as the onset of the afterglow. The single FS component fit is 
consistent with a power law rise index ~3, as expected for an ISM profile and results in fit 
parameters as reported in Tables [2] and [3l This simple model alone cannot explain the steep 
decay of the initial peak. A further emission component at ~2.5 ks superimposed to the 
overall power law rise and decay is required to explain the observed light curve features. Re- 
brightening episodes and variabilities in the optical afte rglow light curve have been observed 
in a number of pr evious burst, e.g. GRB 021004 fe.g. IF ox et al.l |2003[) GRB 030329 (e.g . 
Lipkin et al.ll2004[ ) GRB 050502A JCuidorzi et al.lboosh . GRB 061126 JPerlev et al.lboosl ) 
and GRB 070125 (jUpdike et al.ll2008h . Possible explanations include inhomogeneities iri the 



density profile of the circumburst medium (e.g. IWang fc Loebl 12000 : 
the angular distribution of the outfiow (i.e. the patchy shell mode l, 



or lat e energy injection by refreshed shocks (e.g. iRees fc Meszaroslll998l ). iNakar fc Granot 



jazzati et al. 



2002), in 



Kumar fc Piran 



200C) 



( 120071 ) find that sharp rebrightenings in the optical light curve are very unlikely to be caused 
by density jumps in general and favour the refreshed shock (GRB 030329) or patchy shell 
model (GRB 021004). Points in the afterglow light curve of GRB 070802 which are at- 
tributed to the superimposed component at ~2.5 ks have been excluded from the generic 
light curve fit of the exclusive FS model. 

The peak times of both parametrizations can then be used to constrain the initial 
conditions in the ejecta. Due to the lack of sufficient time resolution between 3 ks and 10 ks, 
the two models cannot be clearly discriminated. 

Within the interpretation of the early brightening of the afterglow of GRB 070802 as the 
RS, a wind-shaped circumburst medium and Poynting fiux dominated models are implicitly 
ruled out. In both cases a RS would not cause a decay and it is very unlikely to be 



imitated by other phenomena flNakar &: Piranll2004l ). 



Although the profile of the light curve of the optical/NIR afterglow shows evidence for a 
reverse shock in an ISM circumburst medium, there is still a problem with its timing relative 
to the burst. The putative RS peak at ~2 ks is highly delayed with respect to the duration 
of the burst (16.4 ± 1.0 s), which is possible for a very thin shell only. As a consequence, 
the initial Lorentz factor of the outfiow, which is estimated to be Fq ~ 40(-E'52/n)^/®, is 



quite small compared to th e expected distribution of the Lorentz f actors above 100 (jPiran 
2000[ ) and previous bursts jMolinari et al.l I2OO7I : IPe'er et all I2OO7I : IPerrero et all 120081 ). A 
possible solution within the context of the RS scenario would be an extremely low density 
environment. 
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Alternatively, the initial incr ease in brightness might be related to the onset of the for- 
ward shock of the aft erglow itself (jPanaitescu fc Kumarll2000l ). as suggested for GRB 060418 
and GRB 060607A (IMolinari et al.l 120071 ) . In this case, a RS component might be hidden 
under the dominating FS emission or occured even before the GROND observations. A pos- 
sible explanation for the superimposed component occurring at th e hght curve peak would 



then be e.g. multiple energy injections by refreshed sh ocks fe.g. Ide Ugarte Postigo et al. 



20051 ) or dense clumps in the circumburst medium (e.g. iGuidorzi et al 



2005). Flares and 



rebrightenings are frequently observed in X-ray afterglow light curves (jO'Brien et al.ll2006l : 
Zhang et al.ll2006l ). The fact that the length of the flare (~1200 s) is compatible with the 
start time of the flare (~1800 s post burst) supports the refreshed shock scenario. 



Using ISari fc PiranI fll999f ). IPanaitescu fc Kumarl fl2000f ) and iMolinari et al.l (120071 ) . we 

160 ( ] , and for the wind shaped case to Fq ~ 



estimate Fq for the ISM case to Fq 

1/4 

80 ■ ■ 



Vo.2n 



E53 



. , both more in line with previous bursts and theoretical predictions. Here E^-^ is 
the isotropic-equivalent energy released in 7-rays in 10^^ erg, 779.2 the 0.2 normalized radiative 
transfer efficiency, n the ISM density in cm~^ and A* the normalized wind density. The 
steep rise of the light curv e with a power law index of 3.56 ± 0.36 favours the homogeneous 
circumburst environment (IPanaitescu fc Kumarl I2OOOI ). 



3.2. The Spectral Energy Distribution 

3.2.1. The intrinsic extinction in previous GRB host galaxies 



The GRB/SN connection hints strongly at the pro genitor of GRBs, which are sup- 
posedly very ma ssive, fast rotating Wolf-Raye t stars (e.g. IWoosley fc BloomI |2006| ) . In this 
collapsar model (IMacFadyen fc Woosleylll999l ) one would expect long GRBs preferentially 
in regions with a high star formation rate (SFR) with significant amount of dust and gas 
in the host along the G RB sight line. However , most previous bursts s how only a mod- 
erate or low reddening (iGalama &: WijersI I2OOII : iKann et al.l l2006l . |2008| ). contrary to the 
observe d high column densi ties of heavy elements and strong depletion of refractory ele- 
ments JSavaglio fc Fallll2004l ) . 



The main feature discriminating between the extinction curves in the Large Magellanic 
Cloud (LMC), Milky Way (MW), Small Magellanic Cloud (SMC) and starburst galaxies is 
the presence and intensity of an absorption feature at 2175 A rest-wavelength. This feature 
is generally associated with the absor ption of graph ite grains, whose abundance and sizes 
changes between the different models (lDraindl2003l ). The feature is most prominent from 
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MW to LMC, whereas it is practically absent in SMC (ICardelli et al.lll989l ) and starburst 
galaxies (jCalzettil l200ll ). SMC and starburst models have a much larger amount of far 
ultraviolet (FUV) extinction. Previous bursts strongly fav oured SMC like dust host galaxies 
( ISchady et al.l 120071 : IStarling et al.ll2007l : iKann et al.ll2008l ). and only for very few of them a 
MW model provides a bette r fit. Additionally , the host extinction shows a trend of lower 
extinction at higher redshift (IKann et al.ll2006l ). However, with increasing redshift, both the 
FUV absorption as well as the 2175 A bump significantly decrease the detection efficiency 
for optical follow-up observations, so the present extinction distribution and its dependence 
on the bursts redshift might be strongly instrumental biased. 



3.2.2. The SED of the afterglow of GRB 070802 



In the present analysis we used MW, LMC and SMC h ke extinction models (jSeatonlll979 



Fitzpatrick fc Massalll986l : lFitzpatricklll986l : |Prevotlll984j ) as templates to fit our multi-band 
data. All of the GROND optical and NIR data were obtained simultaneously at the time 
epoch between 1.5 ks and 3.6 ks post burst. As the generic shape of the early light curve is 
achromatic, we can exclude large effects from an evolving spectrum. The d ata were fit by a 



powe r law and extinction templates from LMC, SMC, and MW in HyperZ (IBolzonella et al. 
2000l ). The amount of dust, power law slope and normalization were free parameters in 



the fit. The redshift w as fixed to the spectroscopic redshift of 2.45 obtained by the VLT 
(IProchaska et al.ll2007l l 



The GROND SED is shown in Fig. [3] and was well fit by the LMC (reduced x^=0.67 
for 4 d.o.f) and the MW (reduced x^=1.96 for 4 d.o.f.) extinction models, while the fit was 
considerably worse for SMC like extinction (reduced x^=4.57 for 4 d.o.f.). The large 
difference between LMC and MW models is mostly due to the g' band magnitude, where 
the error in the GROND data is relatively large. The dust extinction in the GRB host 
in the best fit model is A^°"*=0.9 ± 0.3 mag for MW and A(>°"*=1.8 ± 0.3 mag for LMC 
models. We caution that these values are derived using local extinction curves for a galaxy 
at a redshift of 2.45. The extinction curve and thus the amount of dust reddening could be 
significantly different. Additionally, intervening absorbe rs could contribute to the observed 
dust extinction, which is only resolved by the spectrum (IFynbo et al.l 120071 ). 



We detect a strong absorption feature in the GROND i' band. GROND i' is slightly 
narrower than the SDSS i' band and is located at 7630±537 A. The 2175 A bump at redshift 
2.45 (i.e. at ~7500 A in the observers frame) provides the ideal and obvious candidate for 
this broad absorption feature. The 2175 A bump is the dominating spectral signature of dust 
in the interstellar medium (ISM) in the Milky Way and often attributed to small graphite 
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grains processed by star formation flGordon et al.lll997l : lDraindl2003l : lDuley &: Lazarevll2004j ). 
However, the nature of the bump is not tot ally clear. Its strength varies along different sight 
lines in the Milky Way (ICardelli et al.lll989l ) . Differ ent size distributions or different chemical 
compositions could be the origin of this variation (INatta &: Panagialll984l ). 



At high redshift, the search of the 2175 A feature has always been very difficult. It was 
never clearly detected in single objects, for instance in damped Lyman-a systems (DLAs) 
along quasar (QSO) sight lines. The presence of the bump was excluded from a composite 
spectrum of 37 C a II and Mg II absorbers from Sloan Digital Sky Survey (SDSS) QSOs 
( IWild et al.l l2006l ) . but detected in a combined spec trum of 18 galaxies at l<z<1.5 with 
intermediate-age stellar populations (INoU et al.ll2007l). Further dete ctions outs ide the Local 



Grou p have been suggested e.g. f or GRB 050802 (jSchadv et al.ll2007l) at z=1.71 (IFynbo et al. 



absorption (ISavaglio fc Falll l2004l: 



20051 ) and a galaxy at z=0.83 flMotta et al.l \200'£. Despite the presen ce of strong metal 



Berger et al.l l2006l : iFynbo et al.l l2006l ) and the depletion 



of refractory elements ( Savaglij2006 ) it has not been detected in spectra of previous bursts. 



The afterglow of GRB 070802 shows the clearest presence of the 2175 A dust feature at a 
high redshift so far, where it is detected with GRO ND broad-band photometry and VLT 
spectroscopy as shown in Fig. 5 in lFynbo et al.l (120071 ) and Eliasdottir et al. (in preperation) . 



We can estimate the column de nsity of d i fferen t metals along the sight line for GRB 070802, 
following the approach described in ISavaglio J2006h . The A^""* is directly proportional to the 
dust column density. The dust column density is also proportional to the total metal column 
density. For instance, in a galaxy with MW or LMC like visual extinction, y4y°''*=0.5 mag 
or 0.4 mag is expected for a column density of oxygen of log(iVo)=17.7 or 17.1, respectively. 
Assuming a constant dust-to-metals ratio, our Ay'^^=0.9 mag for MW or 1.8 mag for LMC 
gives log(A'"o) = 18.0 and 17.8. This is the oxygen column density in gas form, i.e. assuming 
that the amount of oxygen locked into dust grains is marginal. In the Galactic ISM, oxy- 
gen i n dust grains is negligible, but up to about 40% in the cool gas (ISavage fc Sembach 
19961 ). If we assume this kind of dust depletion, then the total oxygen column density is 
log(A'"o)=18.4 or 18.2, for MW or LMC, respectively. Using the standard nomenclature for 
hydrogen absorption at X-ray energies and assuming solar metallicity, we conclude that the 
column density of hydrogen log(iVH) is about 21.7 or 21.5, for MW or LMC dust, respectively. 



3.2.3. Combined optical/NIR and X-ray analysis 

The XRT spectrum including all data from 150 s to 405 ks post burst with a total 
exposure time of 56 ks is compatible with an absorbed power law with a photon index r=2.1 
± 0.3 (reduced x^=0.7 for 15 d.o.f.). The foreground hydrogen column density was frozen at 



-li- 



the Galactic value of 3x10^° cm"^ (IDickey fc Lockmanlll990l ). The best fit for A^h at z=2.45 



is 1.60l}jQ X 10^^ cm ^, which corresponds to a total gas plus dust oxygen column density 
of A^o = 7.2tli X IQis cm . Given the large uncertainties, this is consistent with the values 
estimated in the previous section. 

The GROND optical/NIR and XRT X-ray data can be fit together to constrain the 
broad-band afterglow spectrum. We fit the data using a power law with absorptio n in the 



X-ray regime and extinction in the optical/NIR regime using the XSPEC package (lArnaud 



19961 ). Unfortunately, there is no simultaneous early coverage by XRT and GROND, as there 
is a gap in the XRT coverage of GRB 070802 from 911s to 4.2 ks (Fig. [2]). Therefore XRT 
data from 460 s to 5 ks with a total exposure time of 1.3 ks and GROND data from 1.5 ks to 
3.6 ks after trigger were selected for the joint fit. These intervals were chosen because they 
are the closest in time while the source is still bright and do not include the steep decay in 



the early XRT light curve, see Fig. [2] and llmmler et al.l (120071 ). 



The XRT spectrum alone from 460 s to 5 ks post burst is compatible with the total spec- 
trum with a power law photon index of 1.9±0.4, although is no longer well constrained 
due to poorer statistics. 

The combined GROND and XRT SED from 2.2 fim Ks to ~4 keV w as fit by a sin- 



gle po wer l aw with reddening Eb-v a nd A^h Galactic foreground values from ISchlegel et al. 
(119981 ) and lDickey fc Lockmaru (jl990l ) (0.026 and 3xl0^°cm~^, respectively) and by extinc- 



tion in the host at z=2.45. There were insufficient statistics above 4 keV in the XRT band 
to extend the fit to higher energies. The optical/NIR data were modeled using the zdust 
model, where the extinction in the host is obtained using MW, LMC and SMC redde ning 



laws with values of Ry of 3.08, 3.18 and 2.93 respectively, where Ay = Eb-v x Ry (jPei 



19921). As shown in Fig. O the GROND data alone is well described by MW and LMC 
models. The additional absorption in the soft X-ray band was modeled by an absorbing 
column at z=2.45. 

In the broad-band fits we are mainly interested in the power law index between the 
GROND and XRT bands. The combined optical/NIR and X-ray data are well described 
by a single power law with a photon index of 1.91±0.04 (reduced of 1-36 for 12 d.o.f.) 
using the MW extinction model with Eb-v of 0.35±0.04. The A^h at z=2.45 is not well 
constrained by the fit. Similar values of the power law index and Eb-v are achieved using 
the LMC and SMC extinction curves. The results of the zdust and single power law models 
are presented in Table H] with the MW yielding the best fit. 

To test for a cooling break, a broken power law with extinction and soft X-ray absorption 
was also fit to the data. Here the high energy photon index r2 was linked to Fi via F2=Fi+0.5 
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as expected in the fireball model. The value of the break energy was constrained to lie 
between the optical and X-ray bands and the best fit photon indices for the low and high 
energies were ri=1.61±0.05 and r2=2.11 respectively. The reduced for this fit is 1.47 for 
11 d.o.f for the MW extinction curve, 2.40 for 11 d.o.f for LMC and 3.5 for 11 d.o.f for SMC. 
The fit parameters Eb-v and the low energy power law index are presented in Table HI 

The broad-band model without a break provides a better fit to the data, however the 
difference in the reduced parameters of the fits is not conclusive. The single power law 
implies that the cooling break is redwards the GROND bands at this time, which would be 
surprisingly early after the explosion. We caution that the GROND and XRT data are not 
simultaneous. 

We also fit a late time SED using the XRT data 50 ks to 196 ks post burst with the 
GROND data obtained in the JHKs band from 86 ks to 96 ks. The best fit model has a 
power law slope of 2.0±0.2, so the data are still compatible with the early time single power 
law model, but we are not able to distinguish between a single and broken power law fit. We 
do not find evidence for a cooling break between the NIR and X-ray data. 



4. Conclusions 

The optical/NIR afterglow light curve of GRB 070802 can be explained using two mod- 
els. A combined reverse - forward shock model, and a single forward shock model with a 
superimposed emission component at peak brightness. Due to the fact that the afterglow 
peak is heavily delayed compared to the duration of the burst, it is very likely that the in- 
crease in brightness in the early light c urve is related to the onset of the afterglow as proposed 



for GRB 060418 and GRB 060607A flMolinari et all 120071 ). Using an analogue analysis for 

/ \ 1/8 

GRB 070802, we derive an initial bulk Lorentz factor in the iet of around Fn ~ 160 ( — — ) 
' ^ \r10.2n J 

for an ISM environment. The steep rise of the early light curve favours an homogeneous over 
a wind shaped circumburst medium. The ground based observations of the optical/NIR 
afterglow of GRB 070802 were fast enough to detect an early brightening of the afterglow. 
Further rapid follow-up campaigns may establish whether this rapid rise is a generic feature 
of GRB afterglow light curves. 

A broad-band fit of GROND and XRT data is compatible with a single power law 
spectrum with photon index 1.91±0.04, suggesting the cooling break being redwards of the 
GROND bands at the start of the observations at ~1.2 ks post burst. The late time photon 
index from NIR to X-rays is still comparable with the early time power law, indicating no 
time evolution of the spectrum during our observations. 
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The observed SED from g' to the Ks band can be well reproduced with LMC and 
MW extinction models. A broad-band absorption feature in the GROND i' band with a 
central wavelength of 7630 A is unquestionably required to explain the observed SED. The 
redshifted 2175 A feature in the host galaxy of the burst at z=2.45 known from MW and 
LMC extinction models provides the ideal candidate. Depending on the model the best 
fit extinction ranges from y4y'^*=0.9 mag for MW like dust absorption to Ay"^^= 1.8 mag 
for LM C dust. The amount of extinction is significantly large r than estimated for previous 
bursts (IKann et al.ll2006l : ISchady et al.ll2007l : iKann et al.1120081 ). Pre-5'w/it bursts have shown 
a correlation of decreasing extinction with increasing redshift ( IKann et al.ll2008l ). However, 
there might be a strong instrumental bias, as fast and simultaneous optical to NIR follow-up 
observations for a large GRB sample is missing. The GROND instrument, with its unique 
optical and NIR capabilities, is a powerful tool which might remove this bias. GRB 070802 
was the first burst for GROND occurring during nighttime and revealed significant amount 
of dust in its host galaxy, indicating that at least a good fraction of the UVOT dark burst is 
due to intrinsic extinction in the GRB host galaxy. Future observations will help to quantify 
the amount of highly extinguished bursts and may help to solve the mystery of dark bursts. 
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Table 1. Secondary standards in the GRB field in the GROND filter bands used during 

the calibration 



Identifier 


g' 


r' 


i' 


z' 


J 


H 


Ks 




mag 


mag 


mag 


mag 


mag 


mag 


mag 



A 22.42 ± 0.05 20.71 ± 0.04 19.26 ± 0.04 18.50 ± 0.04 16.98 ± 0.06 16.55 ± 0.06 16.09 ± 0.08 

B 16.33 ± 0.07 15.89 ± 0.03 15.78 ± 0.03 15.65 ± 0.03 14.81 ± 0.05 14.59 ± 0.05 14.42 ± 0.06 

C 15.47 ± 0.05 14.98 ± 0.03 14.86 ± 0.03 14.72 ± 0.03 13.87 ± 0.05 13.61 ± 0.05 13.43 ± 0.06 

D 20.07 ± 0.06 18.76 ± 0.04 18.23 ± 0.04 17.92 ± 0.04 16.72 ± 0.07 16.19 ± 0.07 15.74 ± 0.09 

E 15.70 ± 0.05 14.62 ± 0.03 14.17 ± 0.03 13.87 ± 0.03 12.71 ± 0.05 12.61 ± 0.05 11.99 ± 0.06 



Table 2. Parameters of the generic fight curve fit to the r'z'JHK band 



Model 


RS rise index 


RS decline index 


4 

FS rise index 


FS decline index 


ti [s] 


s/ 

Break sharpness 


Combined RS + FS 
Single FS 


3.16 ± 0.33 


-2.66 ± 0.96 


3.12 ± 1.49 
3.56 ± 0.36 


-0.63 ± 0.14 
-0.68 ± 0.04 


2181 ± 107 
1830 ± 51.9 


3.15 ± 2.16 
2.48 ± 0.72 



Table 3. Parameters of the individual fits to the r'z'JHK bands. 



Model 


Parameter 


r' 


z' 


J 


H 


Ks 


Combined RS + FS 
Combined RS + FS 
Combined RS + FS 


F$ [mJ] 
reduced x^/d.o.f. 


14.7 ± 0.7 
4.3 ± 0.7 
0.64/5 


30.9 ± 1.6 
9.7 ± 1.0 
0.98/7 


102 ± 2.5 
33.1 ± 1.6 
1.80/18 


163 ±3.1 
46.5 ± 2.1 
1.10/15 


365 ± 11.0 
135 ± 10.9 
1.53/18 


Single FS 
Single FS 


FO [mJ] 
reduced /d.o.i. 


10.2 ± 0.6 
0.62/3 


21.8 ± 0.9 
0.66/5 


70.6 ± 1.0 
0.98/7 


104 ± 2.0 
0.67/7 


278 ± 7.6 
0.80/7 



Table 4. Broad band spectral fits to the GROND and XRT data using XSPEC. 



Extinction Model Power law E(B-V) Photon Index x^/d-o.f 



MW 


Single 


0.35±0.04 


1.91±0.04 


1.36/12 


MW 


Broken 


0.41±0.04 


1.61±0.05 


1.47/11 


LMC 


Single 


0.35±0.04 


1.92±0.03 


2.13/12 


LMC 


Broken 


0.39±0.04 


1.61±0.05 


2.40/11 


SMC 


Single 




j-yu_o Q4 


3.2/12 


SMC 


Broken 




1.61±0.03 


3.5/11 
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Fig. 1.— NIR afterglow of GRB 070802 in the J band at t=3.4±0.3 ks post burst, including 
the Swift XRT error circle. Only a ~3' x 2.5' segment of the original 10' x 10' image 
is shown. The GROND NIR images have a pixel scale of O'.'6/px each. The image is a 
combination of 48 stacked 10 s exposures and also shows the secondary standards used for 
cahbration, marked as A, B, C, D and E. 
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Fig. 2.— X-ray and optical/NIR light curve of the afterglow GRB 070802. The X-ray 
l ight curve in the upper panel was obtained from the Swift XRT light curve repository 
(lEvans et al.l 120071 ). The afterglow light curve in the GROND bands r', z', J, H and Ks 
is shown in the lower panel. In the g' and i' bands the afterglow was too dim to construct 
a light curve. Also shown are the best fit models in dotted lines for the RS plus FS model 
and in dashed-dotted lines for the FS model. The simple model of an exclusive FS emission 
cannot explain the light curve shape around peak brightness and requires an extra emission 
component superimposed to the overall rise and decay. 
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Fig. 3.— Spectral Energy Distribution of the afterglow of GRB 070802 in the GROND filter 

bands. The data were obtained simultaneously in all colors between To+1.5 ks and To+3.6 ks 
post burst. The shape of the SED can be well reproduced by LMC and MW extinction 
models. A broad extinction feature is detected in the GROND i' band at 7630 A central 
wavelength, which we relate to the 2175 A dust feature in the host at redshift 2.45. Depending 
on the used model, the best fit extinction varies between A^"^* = 0.9 mag (MW) and 1.8 mag 
(LMG). 



